Improved scheme for generation of vibrational trio coherent states of a trapped ion 
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We improve a previously proposed scheme (Phys. Rev. A 66 (2002) 065401) for generating 
vibrational trio coherent states of a trapped ion. The improved version is shown to gain a double 
advantage: (i) it uses only five, instead of eight, lasers and (ii) the generation process can be made 
£0 ' remarkably faster. 
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Following the pair coherent state (PCS) [1], which has proved very important in quantum optics [2], atom-field 
dynamics [3] , quantum-mechanics- versus-local-realism test [4] , continuous- variable quantum nonlocality [5] and quan- 
tum information [6], the so-called trio coherent state (TCS) [7] and their cat-type superpositions [8] have also been 
introduced recently. This novel family of states has been shown inherently nonclassical exhibiting new types of mul- 
timode antibunching, higher-order squeezing and violation of Cauchy-Schwarz inequalities. Although the recognized 
kinds of nonclassical states have been numerous to date (see, e.g. [9]), it is not yet justified that all of them acquire 
£> \ actual applications. It is also not excluded that "really needed" states remain undiscovered or they are among the 
!**""" ■ existing ones with their necessary properties unnoticeable. In that sense, any further detailed study of a known state 
or introduction of a new state would be perhaps equally welcome. As for the TCS, because of their three-mode nature, 
they would play a significant role in phenomena in which the PCS could not. For example, with the TCS one might 
predict that "one-event" continuous-variable local realism violating experiments could be found in a way more or 
less similar to those devised by Greenberger-Horne-Zeilinger [10] or/and by Hardy [11] for observables with discrete 
spectra. However, before any expected experiment can be done, a primary question arises: how can one produce the 
TCS in practice ? This question was answered in [7] for electromagnetic fields interacting with an atom when there is 
f> I strong competition between trio parametric conversion and trio absorption. Here we shall concern with generation of 
(Jj ■ TCS in the center-of-mass motion of a trapped ion. The motivation is that advanced laser-cooling techniques applied 
to single trapped ions allow to cool them down near their zero-point vibrational energy (see, e.g., Ref. [12]). In that 
quantized regime, the ion vibrational modes can be looked upon as ideal bosons which are almost uninfluenced by the 
q-( environment resulting in negligible decoherences. This fact greatly favors potential implementations such as designing 
!• ■ quantum logic gates in a quantum computer [13] (in particular, the controlled-NOT gate has been realized in two 
40 Ca + ions in a Paul trap individually radiated by focused lasers [14]). Laser-driven trapped ions placed inside optical 
cavities can also entangle the cavity mode simultaneously with the internal and external degrees of freedom of the 
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ion [15]. Previously, the vibrational TCS in a 3D trap was shown to be generated by a scheme using eight lasers [16] 
(compare with schemes for the vibrational PCS in a 2D trap [17]). In this Letter we shall improve the scheme proposed 
in [16] by reducing the number of used lasers from eight to five. Moreover, we shall carry out detailed simulation to 
reveal that the time needed to generate the vibrational TCS can be made shorter in the improved version than in the 
original one [16]. 

The vibrational TCS, denoted by |£, p, q) with £ = r exp(iip) and p, q integers, is a three-mode entangled continuous 
superposition of coherent states of the form 
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where N(p,q,r) defined by 



r 2n 



N(p,q,r)- 2 = V (2) 
^- / {n + q)\(n+p)\n\ 

is the normalization coefficient and |--)x(j/,z) stand for coherent states in the vibrational mode along the x (y,z) 
axis with the annihilation boson operators a x ( y ^ z y Equivalently, the state \£,p, q) can also be represented in terms of 
correlated boson number trios as 

oc 

\LP,q) = ^2C n (Z,p,q)\n + q) x \n + p) y \n) z (3) 

n=0 

where 

Cn K, P, q) = . — (4) 

y/(n + q)\(n+p)\n\ 

and \---) x t yz \ stand for Fock states in the x (y,z) axis mode. Both the representations (1) and (3) come from the 
fact that the vibrational TCS is defined as the joint eigenstate of the operators a x a y a z , P = a+a y — afa z and 
Q = a+a x - a+a z , i.e. 

a x a y a z \€,p,q} =€\£,p,q) , (5) 
P\S,P,q)=p\Z,P,q), (6) 

Q \€,p,q) = q \€,p,q) ■ (7) 

Generally, the system of a two-level (\g) and \e) : ground and excited state) ion localized in a small region by a 3D 
isotropic harmonic potential and radiated by L external lasers is described by the Hamiltonian (U = 1) 

H = H + Hi n t, (8) 
H = Ia<j z + V v(a}a j + -), (9) 



j=x,y,z 



H in t - I 0i ex P[*(^* + <f>l)]9l<r- + H c -] ( 10 ) 

with A the energy gap between the ion's two levels, er_ = \g) (e| = a z = |e) (e| — |^} (g\ , v the energy of a 
quantum of the vibrational motion, f2; the Rabi frequencies, u>i ((pi) the laser frequencies (phases) and gi the laser 
spatial profile. For travelling waves gi — exp[— ik;R;] with R; the position operators along the laser propagation 
directions determined by the wave vectors k;. Instead of 8 lasers as proposed in [16], here we use only 5 lasers. The 
first four lasers are directed as sketched in Fig. 1, each of which is detuned to the third lower sideband of the ion 
vibrational motion, i.e. w; = i,2.3,4 = A — 3v. As for the fifth laser, its direction is unimportant but its frequency 
must be in resonant with the ion transition, i.e. w; = 5 = A. Assuming an equal wave number for all the lasers we 
are left with a unique Lamb-Dicke parameter rj which is a measure of the ion's localization region compared with the 
laser wavelength. If v is much larger than any other characteristic frequencies (resolved sideband limit) and r\ <C 1 
(Lamb-Dicke limit), we have, in the interaction picture and in leading order in rj, the interaction Hamiltonian of the 
form [16] 



H 



int 



i=i 



(T++H.C. (11) 
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where the annihilation operators of the vibration quanta in the direction R/ are denoted by Ai which are related to 

Ai = a x +a y + a z , A 2 — a x — a y + a z , (12) 

A3 = a x +a y -a z , A 4 = a x — a y -a z . (13) 

If we now adjust the laser intensity and phase to meet the conditions f2i=i,2,3,4 = ^ and ^=1,3 = $1=2.4 + t = <f>, then 
the Hamiltonian (11) simplifies to 

Hint = C (a x a y a z - £) cr+ + H.c. 
with £ and £ determined by the system parameters in a controllable manner as 

( = -6n 1 3 nexp[-{i$ + V 2 /2)} 7 (14) 



e = -en^«p[<w-^)]- (15) 



Because the damping of vibrational quanta is negligible, the main decay process occurs via ionic spontaneous emission 
with rate 7. Thus, in the Lamb-Dicke limit, the system density operator p obeys the master equation [16] 

dp 'y 

— = -i[Hi„ t ,p] - - {<?+o-P + po+o- - 2a-po + ) . (16) 

In the long-time limit, Eq. (16) has a "dark" steady solution p s = \g) \^f) (W| (g\ with \^>) responsible for the steady 
state of the ion vibration. It is easy to verify that obeys the equation 

a x a y a z \*)=Z\*). (17) 

Further, if the ion is initially prepared in a Fock state \^o) = 1 9 + k) x \p + k) \k) z with k a non- negative integer, 
then, because p and q are conserved in the problem under treatment, the state f^) satisfies also the two equations 

P|*)=p|*>, (18) 

Q|*) = g|*>. (19) 

Comparing (17) - (19) with (5) - (7) yields \^) = \£,p, q) which is the vibrational TCS we want to generate. 

To follow the transient process and to assess generation time (i.e. the time needed for the system to reach the 
steady regime) we simulate Eq. (16) by the Monte Carlo Wave- Function approach [18]. The system wave function at 
time t is looked for in the form 

$(t) - Yl [ G ""«W l3) + E 'mn(t) |e>] \l) x \m) y \n) z . (20) 

For the purpose of numerical calculation, the Fock-state basis should be confined to l 7 m 7 n < N max with iV max a 
cutoff to be chosen such that its increasing does not change the result within a pre-set high enough accuracy. In 
our simulations we have found iV max = 10 sufficient for most of the numerical calculation: just an error of less 
than 10 -6 arises by increasing the cutoff by one. For convenience, we use the dimcnsionlcss time r determined by 
t = -ft. For the initial condition $(0) = |e) \q + k) x \p + k) y \k) z with k = 4, the time evolution of the probability 
P(l, m, n) — |G;„ m | 2 + \Ei mn \ 2 of finding the ion in a Fock state \l) x \m) y \n) z is plotted in Fig. 2a. As is visual from 
that figure, at the beginning P(l,m,n) oscillate but at large r they approach \C n (£,p, q)\ 2 (see Eq. (4)) precisely, 
implying generation of the vibrational TCS. The initial variation of P(l,m,n) depends on the value of k but their 
long-time behavior is /c-independent. This means that the TCS is eventually generated independently of the choice of k 
which is best taken as k — because the state e) \q) x \p) |0) z is easiest to prepare practically. The transition towards 
the steady regime can be watched as well in Fig. 2b showing how the inversion (<j z ) = J2i m n (\^imn\ 2 — |G; m „| 2 ) 
evolves from 1 to —1. The purity of the generated state is examined by the fidelity F(t) = | ($(t)| £,p,q) | 2 whose 
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change in time is depicted in Fig. 3. Because of the truncation in the Fock basis we determine the vibrational TCS 
generation time r s by the constraint 1 — F(t s ) = 10~ 5 and display the r s as a function of a = C/7 in Fig. 4a. The 
decreasing of t s with increasing a plays a role as will be seen. 

Finally, we discuss a possibility to make the generation time in the present improved scheme (L = 5 scheme) shorter 
than that in the previous scheme (L = 8 scheme) [16]. Suppose we wish to produce a state q) with fixed £, p and 
q. The desired values of p and q are easily to manage by preparing an initial state in the form |e) \q + n) x \p + n) y \n) z 
with n a non-negative integer. To tailor £ we use Eq. (10) of Ref. [16] and Eq. (15) of this work. The same value 
of £ in the L = 8 scheme can be obtained in the L = 5 scheme in either of the two following ways: (i) keeping the 
frequency of the 5th laser in the L = 5 scheme equal to that of the 8th laser in the L = 8 scheme, but decreasing 
the frequencies of the 1st to 4th laser in the L = 5 scheme 6 times smaller than those of the 1st to 7th laser in the 
L = 8 scheme or (ii) increasing the frequency of the 5th laser in the L = 5 scheme 6 times larger than that of the 
8th laser in the L = 8 scheme, but keeping the frequencies of the 1st to 4th laser in the L = 5 scheme equal to those 
of the 1st to 7th laser in the L = 8 scheme. In case (i) the value of ( in the L = 5 scheme (see Eq. (14) of this 
work) remains unchanged in comparison with that in the L — 8 scheme (see the equation before Eq. (10) in Ref. 
[16]) and, hence, both schemes undergo the same dynamics. Nevertheless, if case (ii) is chosen, the value of ( in the 
L = 5 scheme becomes 6 time larger than that in the L = 8 scheme. Since t s , for the same value of 7, decreases with 
increasing ( (see Fig. 4a), a net advantage in the generation time appears in the improved scheme. Figure 4b shows 
that advantage from an angle of the ^-dependence, for more clarity. A noteworthy merit of the L = 5 scheme is that, 
for the whole range of £, the desired TCS is generated in a remarkably shorter time compared with the L = 8 scheme. 

In summary, we have improved a recently proposed scheme for generating trio coherent states in the vibrational 
motion of an ion trapped in three dimensions. The trio coherent state generated by this scheme is stable because it 
appears in a steady regime in which the ion has fully relaxed to its ground "dark" state (Fig. 2b) which is identified 
by the quenched fluorescence from the ion. Besides reduction in the number of driving lasers to be used, the improved 
version can also provide a much shorter generation time. The possibility to shorten the generation time t s relies 
upon its decrease with increasing a. In general, such a decreasing behavior does not exist for the whole range of a. 
In fact, t s would turn out to increase with a, say, when a is very large. For instance, in the extreme limit of a — > 00 
(corresponding to negligible spontaneous emission, 7 — > 0), periodic oscillations or phenomena like collapses and 
revivals, undamped self-pulsing or even chaos would occur. However, in problems like the one under consideration, 
where relaxation to a steady state within a finite time is important, 7 is essential and cannot be treated as arbitrarily 
small and, thus, arbitrarily large values of a are irrelevant. To specify the range of a for which our simulations have 
been carried out, we use 77 as a measure of smallness since 77 -C 1 in the Lamb-Dickc limit we are working in. As 
explained, because the ionic decay cannot be neglected, let us be interested in situations such that 7 satisfies the 
condition 7/(60) 3> 0(rf). Then, from Eq. (14), we get a — C/7 ~ 0(6ili] 3 /-/) which is much less than unity. For 
smaller values of 7 (i.e. larger a), as mentioned above, abundant dynamics associated with other different timescales 
will take place which should be dealt with separately (e.g., by means of multiple-timescale analyses). 
Acknowledgments 

The authors thank the KIAS Quantum Information Group for useful discussions. H.S.Y. is supported by a KAIST 
Grant for Basic Science, B.A.N, by a KIAS Research Fund (No. 02-0149-001) and J.K. by a Korea Research Foundation 
Grant (KRF-2002-070-C00029). 



[1] G. S. Agarwal, Phys. Rev. Lett. 57 (1986) 827; 

G. S. Agarwal, J. Opt. Soc. Am. B 5 (1988) 1940. 
[2] Ts. Gantsog and R. Tanas, Opt. Commun. 82 (1991) 145; 

B. A. Bambah and G. S. Agarwal, Phys. Rev. A 51 (1995) 1698; 

A. Gilchrist and W. J. Munro, J. Opt. B: Quantum Semiclass. Opt. 2 (2000) 47; 

Arvind, Phys. Lett. A 299 (2002) 461. 
[3] A. Joshi and R. R. Puri, Phys. Rev. A 42 (1990) 4336; 

S. C. Gou, Phys. Rev. A 48 (1993) 3233. 
[4] K. Tara and G. S. Agarwal, Phys. Rev. A 50 (1994) 2870; 

A. Gilchrist, P. Deuar and M. D. Reid, Phys. Rev. Lett. 80 (1998) 3169; 

W. J. Munro, Phys. Rev. A 59 (1999) 4197; 

A. Gilchrist, P. Deuar and M. D. Reid, Phys. Rev. A 60 (1999) 4259. 



4 



[5] S. Mancini, V. I. Manko, E. V. Shchukin and P. Tombesi, e-print quant-ph/0302089. 

[6] S. Mancini and P. Tombesi, Quant. Inf. and Comp. 3 (2003) 106. 

[7] B. A. Nguyen and M. D. Truong, J. Opt. B: Quantum Semiclass. Opt. 4 (2002) 80; 

B. A. Nguyen, J. Opt. B: Quantum Semiclass. Opt. 4 (2002) 222. 
[8] B. A. Nguyen and M. D. Truong, J. Opt. B: Quantum Semiclass. Opt. 4 (2002) 289; 

B. A. Nguyen, Phys. Lett. A 311 (2003) in press. 
[9] V. V. Dodonov, J. Opt. B: Quantum Semiclass. Opt. 4 (2002) Rl. 
[10] D. M. Greenberger, M. A. Home and A. Zeilinger, in Bell's Theorem, Quantum Theory and Conceptions of the Universe, 

ed. M. Kafatos (Kluwer Academic, Dordrecht, 1989); 

N. D. Mermin, Am. J. Phys. 58 (1990) 731; 

B. Yurke and D. Stoler, Phys. Rev. Lett. 68 (1992) 1251. 
[11] L. Hardy, Phys. Rev. Lett. 68 (1992) 2981; 

L. Hardy, Phys. Lett. A 167 (1992) 17. 
[12] C. Monroe et al, Phys. Rev. Lett. 75 (1995) 4011. 
[13] J. I. Cirac and P. Zoller, Phys. Rev. Lett. 74 (1995) 4091; 

C. Monroe et al., Phys. Rev. A 55 (1997) R2489; 

F. Poyatos, J. I. Cirac and P. Zoller, Phys. Rev. Lett. 81 (1998) 1322; 

J. A. Sorensen and K. Molmer, Phys. Rev. A 62 (2000) 022311; 

A. M. Chids and I. L. Chuang, Phys. Rev. A 63 (2000) 012306. 
[14] F. Schmidt-Kaler et at, Nature 422 (2003) 408. 
[15] A. B. Mundt et al, Appl. Phys. B 76 (2003) 117; 

S. S. Sharma, Phys. Lett. A 311 (2003) 111. 
[16] B. A. Nguyen and M. D. Truong, Phys. Rev. A 66 (2002) 065401. 
[17] S. C. Gou, J. Steinbach and P. L. Knight, Phys. Rev. A 54 (1996) R1014; 

W. Kaige et at, Phys. Rev. A 63 (2001) 043419. 
[18] J. Dalibard, Y. Castin and K. Molmer, Phys. Rev. Lett. 68 (1992) 580. 



5 



Figure captions 

Fig. 1: The configuration of used lasers : the 1st (2nd, 3rd and 4th) laser propagates along the direction connecting 
the coordinate origin "o" with the point with {x, y, z} = {1, 1, 1} ({1, —1, 1}, {1, 1, —1} and {1, —1, —1}), while 
the 5th laser can be oriented arbitrarily (not shown). 

Fig. 2: a) Left: The probability P(l,m 7 n) versus scaled time t = jt for the initial condition $(0) = 
|e) \q + A) x \p + A) y |4) z . The parameters used are a = 0.02, £ = 2.0, p = 3 and q = 2. The values of I, m, n 
are indicated by (l,m,n). All the P(l,m,n) with m ^= n + p and I =/= n + q are identically zeros, whereas all 
the P(q + n,p + n,n) with n > 3 are negligible in the long-time limit. Right: Distribution of |C„(£,p, q)\ 2 (see 
Eqs. (3) and (4)) for the same values of £, p and q as above. The values of n are indicated by (n). All the 
\C n (l;,p,q)\ 2 with n > 3 are negligibly small, b) Time evolution of the ionic inversion (a z ) for the same initial 
condition and parameters as in a). The inset enlarges the variation for the interval r G [0, 20]. 

Fig. 3: Fidelity F as a function of r for £ = 2.0, p = 3, q = 2 and a = 0.01, 0.02, 0.03 (curves indicated by "a", "b", 
"c", respectively). 

Fig. 4: a) The generation time r s as a function of a for £ = 2.0, p = 3 and q = 2. b) As in a) but in logarithmic scale 
and as a function of £ in the L = 8 and L — 5 schemes. Here we use a = 0.02 in the L = 8 scheme which 
corresponds to a = 0.12 in the L = 5 scheme. 
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